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ABSTRACT 

By examining the diffusion of young white dwarfs through the core of the globular cluster 47 Tu- 
canae, we estimate the time when the progenitor star lost the bulk of its mass to become a white dwarf. 
According to stellar evolution models of the white-dwarf progenitors in 47 Tucanae, we find this epoch 
to coincide approximately with the star ascending the asymptotic giant branch (3.0 ± 8.1 Myr before 
the tip of the AGB) and more than ninety million years after the helium flash (with ninety-percent 
confidence). From the diffusion of the young white dwarfs we can exclude the hypothesis that the 
bulk of the mass loss occurs on the red-giant branch at the four-sigma level. Furthermore, we find 
that the radial distribution of horizontal branch stars is consistent with that of the red-giant stars 
and upper-main-sequence stars and inconsistent with the loss of more than 0.2 solar masses on the 
red-giant branch at the six-sigma level. 

Subject headings: globular clusters: individual (47 Tuc) — stars: Population H, Hertzsprung-Russell 
and C-M diagrams, kinematics and dynamics 


1. INTRODUCTION 

When and how a star like the Sun loses its mass to 
become a white dwarf star is a key open question of 
stellar evolution. Does the bulk of the mass loss occur 
when the star is a red giant or when the star reaches the 
asymptotic giant branch? For forty years, in stellar evo¬ 
lution models mass loss on the red-gian t branch ( RGB 
has been commonly described with the Reimersj (1975| 
formula, an empirical scaling relation that involves ba- 
sic stellar parameters and an adjustable efficiency pa¬ 
rameter, 77 rgb- This parameter is usually constrained 
from the requirement of reproducing the extended blue 
horizontal branches (HB) in the Her tzsprung-Russell dia¬ 
gram s of Galactic globular clusters (Renzini & Fusi Peed 
19881. Typical values are ryRCB — 0.35 — 0.45, which im¬ 
plies that ~ 0.20 — 0.25 Mq is the mass that should be 
lost on the red-giant branch by a star with initial mass 
0.85 — 0.90M0, characteristic of the turn-offs in Galac¬ 
tic globular clusters. An additional mass loss event of 
smaller size, ~ 0.10 — 0.15Mq, should take place later, 
during the evolu tion on the asympto t ic gia nt branch 
(AGB). Recently McDonald & Zijlstra (2015) have also 
argued for values of ?7 gb = 0.477± 0.070^0^^ for Galac¬ 
tic clusters and ?7 gb = 0.452 for 47 Tuc in particular. 

On the other hand this classical paradigm has been 
recently challenged from various perspectives. Perhaps, 
one of the most intriguing findings is that the extended 
blue horizontal branch may be well explained by the very 
high helium abundance associated with the bluest main 
sequence of the multiple populations widely present in 
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Galactic globu lar clusters (e.g., Lee et al. 2005 D’Antona 
&: Galoi|2008). In general, stellar models indicate that at 


the same age, higher helium corresponds to lower turn¬ 


off mass (Bertelli et al. 


2008). Therefore, for the same 


mass ejected on the red-giarit branch, a higher helium 
abundance favors the development of more extended hor¬ 
izontal branches. 

At the same time there are hints that the mass loss 
rates predicted for red-giant-branch stars with the clas¬ 
sical calibrations may be significantly overestimated. Us¬ 
ing chromospheric model calculations of the Ha line for 
a sa mple of red-giant-bran ch stars in a few globular clus¬ 
ters iMiszSrosTtaL] (p009|) pointed out that the resulting 


mass loss rates are about one order of magnitude lower 
than obtained with the Reimers law or inferred from the 


infrare d excess of similar s tars by Origlia et al. (2007 


see also|Origlia et ari(|2014| for recent results). Such dis- 

Tb 


crepancy a ppears to be overcome according to|Groenewe-| 


gen (2012), who found agreement between the mass loss 
rates derived from fitting the spectral energy distribu¬ 
tions of red-giant-branch stars with infrared excess, and 


the chrom ospheric estimates. More recently Groenewe- 
gen (2014) has performed the first detection of rotational 


GO line emission in a nearby red giant branch with a lu¬ 
minosity of ~ 1300 7^0 and an estimated mass-loss rate 


as low as a few 10 ^ Mq yr ^. Interestingly Miglio et al. 


(2012) argued from Kepler asteroseismic measurements 
of the stars in the very metal-rich open cluster NGG 6971 
that low values of ij on the red-giant branch (?7rgb) are 
needed to account for the mass loss between the red giant 
and red clump phases of stars in this metal-rich cluster 

(0-1 ^ Vrgb ^ 0-3). 

On the other hand, mass loss on the asymptotic giant 
branch of old stellar populations may have been under¬ 
estimated at lower metallicities. Recent stellar popula¬ 
tion synthesis studies have shown that to reproduce the 
star counts and luminosity functions of metal-poor low- 
mass thermally-pulsing asymptotic-giant branch (TP- 
AGB) stars in a sample of nearby galaxies one has to in- 
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voke a more efficient mass loss than the classical Reim ers 
recipe ([Girardi et al. |2010 Rosenfreld et al. 2014 


see 


Sect. I'his also yields good agreement with the low- 
mass end of the initial-fin al mass relation, as probed with 
the white dwarfs in M4 (Kalirai et al.||2009|). 

In t his study we build upon the results of [Heyl et al.| 
(20151 to demonstrate that the bulk of the mass loss from 


stars m 47 Tucanae must happen on the asymptotic gi¬ 
ant branch. In the core of the globular cluster 47 Tu¬ 
canae the timescale for dynamical relaxation through 
two-body interactions is similar to the stellar evolution 
timescale for a star to live as a horizontal-branch star, 
rise up the asymptotic branch and b ecome a white-dwarf 
star ( Harris|[l996 Heyl et ^|2015 |. The core of 47 Tuc 


has been the focus of numerous previous investigations 


(e.K IMcLaughlin et al. 

2006 Knigge et al. 

2008 

Berg- 

busch &: Stetson 2009 

Milone et al. 2012a|, but Heyl 

et al. (|2015|) is the hrst paper that combines the near 


ultraviolet filters of the Hubble Space Telescope (HST) 
with a mosaic that covers the entire core of the cluster. 
In these filters the young white dwarfs, the giant stars, 
the horizontal-branch stars, the blue stragglers and the 
upper main sequ( 
as shown in Fig. 
lar resolution in t 
WFC3/UVIS cam 
and incompleteness in this crowded field. From theoret¬ 
ical arguments it has long been argued that two-body 
interactions will sort the stars in a globular cluster by 
mass with the more mas sive stars lyin g closer to the 


To construct the stellar evolution models here, we used 
both M ESA (Modules for E xperiments in Stellar Astro¬ 
physics; Paxton et al.||2011 ) and a combin ation of PAR¬ 
SEC (for the evolutioii before the TP-AGB ;||Bressan et al. ' 
2012 1 and COLIBRI (for the TP-AGB; Mango et al. 
20131 to perform simulations of stellar evolution starting 
with a pre-main-sequence model of 0.9 solar masses and 
a metallicity of Z = 4x 10~ ^ and Y = 0.256, appropriat e 
for the cluster 47 Tucanae (jBergbusch & Stetson||2009|), 
assuming [Fe/H] = —0.83 and [a/h'ej = 0.3. We note 
that adopting [Fe/H] = —0.76 and individual elemental 
abundances measure d in 47 Tuc stars (as summarized in 
Milone et al. 2012a) leads to a somewhat larger metal- 
licity (Z ~ 0.0055), but the general trends, discussed 
below, do not change significantly. 

For the MESA models, we used SVN revision 5456 
and started with the model lM_pre_ms_to_wd in the 
test suite. We changed the parameters initialunass 
and initial_z of the star and adjusted the parameter 
log_L_lower_limit to —6 so the simulations run well 
into the white dwarf cooling regime. We also reduced 
the values of the wind rj on the RGB and AGB to 0.46 


Furthermore, the exquisite angu- 

dwarf (jRenzini & Fusi Pecci 

1988 

Moehler et al. 

2004 

he near ultraviolet of the new Hubble 

Kalirai et al. 

2UU9) trom the U.9 solar mass progenitor. 


of the tip of the TP-AGB (Lebzelter & Wood 2005) which 
is sensitive to the mass of the resulting W U a nd with the 


center of the cluster (e.g. Spitzer 1987), and this mass 
segregation has been q u aiitified in various clusters (e.g 


2015 

the 


obser ved cooling curve of the white dwarfs (Heyl et al. 
^ . For the mas s loss on the red-giant branch we used 
R( 


.eimers (1975) value 


Goldsbury et al. 2013). Heyl et al. (2015) for the first 
time caught this process of mass segregation in action 
and determined the timescale for the sorting of stars by 
mass, the relaxation time, to be about 30 Myr. They 
have outlined in detail how to measure the completeness 
rate and model the diffusion of stars in the cluster using 
the young white dwarfs. 

Here we build upon this diffusion model by comparing 
the radial distribution of the white dwarfs with that of 
the upper-main-sequence and red-giant stars. The key 
observation that one can draw from Fig. is that the 
distribution of the bright white dwarfs, "vraose median 
age along the white-dwarf cooling track is 6 Myr, more 
closely resembles that of the upper-main-sequence stars 
than that of the fainter white dwarfs of about 100 Myr. If 
the progenitors of the white dwarfs lost their mass more 
than 30 Myr before the birth of the white dwarfs, the 
white dwarfs would have already been sorted by mass, so 
their radial distribution would not look so similar to that 
of the upper-main sequence stars. Furthermore, because 
the horizontal branch in 47 Tucanae is thought to last for 
about 80 Myr, a few relaxation times, their radial distri¬ 
bution will also reflect the mass of the horizontal branch 
stars. The radial distribution of the horizontal branch 
stars is very similar to that of the upper-main-sequence 
stars, their progenitors. We will confront these obser¬ 
vations with expectations from stellar evolution models 
and better quant ify it with the diffusion models from 


Mii = 4x IQ-^^r] 


L R Mq 
Lq Rq m 


[MqY- 




( 1 ) 


and on the asymptotic-giant branch we use the Blocker 
(|1995| formula 


Ms=4.83xl0-®(^) (^) 


( 2 ) 


Heyl et al. (2015). 


2. STELLAR EVOLUTION MODELS 


These parameters yield a model where the star loses 
about 0.2 solar masses as a red giant and 0.17 solar 
masses as an asymptotic giant star. As other options 
we also used a value of rj on the red-giant br anch of 0.1, 
0.2 and 0.3 in the range of Miglio et al. (20121) and higher 
values on the AGB (77 agb) as outlined in 'Tabled These 
also yielded a 0.53 solar mass white dwarf but wrHi much 
less mass loss on the red-giant branc h with values of n i n 
better acj;»rdance with the results of Miglio et al. (2012). 

Table m summarizes the results of the various wind 
models using MESA. Essentially the two wind param¬ 
eters 77RGB and tjagb can be tuned to change the ratio 
of mass loss on the two giant branches without chang¬ 
ing the initial or final mass of the star (here 0.9 and 
O. 53 M 0 ), the age of the star where it becomes a white 
dwarf at the tip of the AGB (two) — this should equal 
the age of the globular cluster today. We examine two 
other timescales. The first is the time interval between 
the tip of the red-giant branch and the tip of the asymp¬ 
totic giant branch, Atx 2 Ti ^cid the second is the time 
between the helium flash and when the central helium 
abundance drops below 10“®, AtcHB, the epoch of core 
helium burning. During the core helium burning stage 
of the star, the luminosity remains nearly constant for 
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Fig. 1.— Left: Color-magnitude diagram of the core of 47 Tucanae in the WFC3 filters F225W and F336W. We have imposed a mild cut 
on magnitude error to remove strong saturated stars from the sample. The inset gives the number of stars in each class before correcting 
for completeness. The classes are: BSS, blue straggler stars; UMS, upper-main-sequence stars of about O.QMq; RGB, red-giant-branch 
stars; HB, horizontal-branch stars; BWD, bright white dwarfs; FWD, faint white dwarfs; and LMS, lower-main-sequence stars of about 
0.65Mq. We determine locations of these classes on the color-magnitude diagram from the PARSEC isochro nes and the MESA evoluti onary 
models. The magenta curves trace evolutionary tracks calculated with MESA and the atmosphere models of|Castelli &: Kurucz|(|2004||. The 
magenta curves from bottom to top are for initial masses of 1.1, 1.4 and 1.8 solar mas ses. We have assu med tliat the post-mam sequence 
of the bli^-straggler stars follows a similar track to normal stars of a similar mass (|Sills et al.||2009 . The green curv e is a PARSEC 
isochron^^ for an age of 11 Gyr and the metallacity of 47 Tucanae using the bolometnc corrections ot Chen et al.| (|2014||. An additional 
0.4 and (k3 magnitudes of extinction were added in F225W and F336W respectively to fit the isochrone to tlie data. The dashed line 
traces F3S6W = 15.25, the approximate magnitude where the images are saturated, so model curves to the right of this line will have 
observed colors just to the left of the line. Right: The completeness-corrected radial distribution of stars in the various regions highlighted 
in the color-magnitude diagram. The radial distribution of the upper-main sequence stars (UMS) and the red giants (RGB) are nearly 
indistinguishable. Here the inset gives the size of each sample after correcting for completeness. 

^http: //stev.oapd. inaf . it/cgi-biii/cmd_2.7 
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TABLE 1 

Wind Parameters for the Various MESA Models. 


Vrgb 

»?AGB 

A7hb 

AtwD 

[Me] 

AtT2T 

[Myr] 

^iCHB 

[Myr] 

IWD 

[Gyr] 

0.1 

0.7 

0.86 

0.54 

108 

89.2 

10.8 

0.2 

0.6 

0.82 

0.54 

109 

89.7 

10.8 

0.3 

0.5 

0.78 

0.53 

110 

90.4 

10.8 

0.46 

0.46 

0.70 

0.53 

111 

91.0 

10.9 


about 80 Myr, somewhat less than these two timescales, 
according to the PARSEC evolutionary tracks. This is 
the expected time that stars will linger within the region 
of the CMD denoted as the horizontal branch in Fig. 
When the two wind parameters are equal, the mass loss 
on the two branches also ends up being about equal and 
using the range of parameters outlined by |Miglio et al.| 
(20121 one can have as little as one-ninth of the mass loss 
on the red-giant branch, leaving nearly ninety percent of 
the mass loss to occur on the asymptotic giant branch. 

The procedure for the PARSEC-COLIBRI models was 
similar, but for the mass-loss descriptions. A first set of 
models was computed adoptin g the Reimers law on the 
RGB, in combination with the Schroder fc Cunt^ ( 2005|) 
formula for the TP-AGB, as modified by|Rosenhela et al. 

((Mil): 


AfmSC = 10 -7^ 


L R Mq / Teff 


LqRq M V4000K 


8.9 


1-L 


g0 \ 
43005 ) 


(3) 


[Mqyv 


where g = GM/R^ is the surface gravity, and 5mSC 
is a free efficiency parameter. Again for each value of 
gRGB = 0.1, 0.2,0.3, 0.4, the parameter 5mSC of the AGB 
mass loss was tuned so as to obtain a final G-0 core with 
mass of ~ 0.53 — 0.54M©. We also found that the tip of 
the TP-AGB in PARSEC-GOLIBRI models was consis¬ 
tent wit h the observed My at th e tip of the TP-AGB in 
47 Tuc ( Lebzelter fc W^ d| |2005 1. 

We recall that Eq. nr was proposed by |RosenfieTd| 
et al.| (|2014D with the purpose of reproduciiig the star 
counts and luminosity functions of TP-AGB stars de¬ 
tected in a sample of nearby metal-poor galaxies from 
the AGS Nearby G alaxy Survey Treasury (ANGST; [Dab 
canton et al.||2009|) , which do not present signs of recent 
star formation. The working scenario is that at lower 
luminosities on the AGB winds are not driven by dust 
opacity, rather they should be linked to the mec hanical 
flux produced in h ighly turbulent chromospheres (Cran- 
mer & Saar 2011|). Adopting 5 rgb = 0.2 for RGB stars, 
the AJNGS'l data were best reproduced assuming that on 
the TP-AGB low-metallicity low-mass stars experience 
more mass loss than predicted by the classical Reimers 
law (with 5 rgb = 0.3 — 0.4). Good agreement with 
observations was obtained by settin g the efficiency pa¬ 
ramet er 77mSC = 0.4 in the modified Schroder fc Guntz 
(20051 formula. The calibration also sets constraints on 
the i'P-AGB lifetimes for [Fe/H] < —0.9. They should 
be ~ 0.5 Myr for lower mass stars (M < IMq), corre¬ 
sponding to final masses of 0.52 — 0.54M©. 

A second test calculation was performed with the pre¬ 
scriptions recently proposed by 
describe mass loss on both the 


Origlia et al. (2014) to 
RGB and the ITOD of 


Galactic globular clusters. This formulation was derived 
relating the mid-IR excess exhibited by a fraction of RGB 
stars to the mass-loss rate, through a scaling relation that 
involves a few parameters for the dust properties. We set 
the parameters equal to the reference values suggested by 
the authors for 47 Tucanae, namely: expansion velocity 
Uexp = 10 km s“^, gas-to-dust ratio S = 200, and grain 
density p = 3 g cm“^. Another quantity to be specified 
is the fraction /on of dusty RGB stars, which may vary 
with the bolometri c magnitude Mroi. Its value was de¬ 
rived from Fig. 4 of Origlia et al. (2014) paper. Assuming 
[Fe/H] = —0.7 for 47 'i'ucanae, we got /on = 0.098 for 
— 1.5 < Mboi < —0.6 and /on = 0.222 for Mboi < —1.5 
on the RGB; /on = 0.292 for Mboi < —1.5 on the AGB. 
No mass loss was considered for Mboi > —0.6. 

A third test calculation was carried out using another 
semi-empirical mass loss relation based on the measured 
infrared excess. With the aid of dust radiative transfer 
calculations that best fit the observed spectral energy 
distributions of a sample of field RG B stars, with accu¬ 
rate parallaxes, Groenewegen (2012) derived the formula 


[Meyr-^j . (4) 


We note that there is no adjustable param eter here, and 
that th e reference dust parameters used by |Groenewegen| 
(2012) were 5 = 200, and Ugxp = 10 km s which are 
exactl y the same as the ones adopted in the |Origlia et al.| 
(2014) formulation. Instead, as we will see below, the 
predictions in terms of RGB mass loss are significantly 
different! Later on the AGB mass loss was described 


with the modified Schroder & Guntz (20051 formula given 
by Eq. (|^. An efficiency parameter PmSC = 0.8 was 
found to oe the suitable choice to obtain a final mass of 
~ 0.54Mo. 

Table [^outlines the results of the PARSEC-COLIBRI 
models. The quantitative trends are quite similar to 
those obtained with the MESA code. Small differences in 
the mass on the horizontal branch for the same prgb can 
be easily explained by the details of the mo del. Concern¬ 
ing the model with the Origlia et al. (2014) mass loss, we 
just note that it yields a hnai mass of 0.57 Mq, that is 
larg er than our re f erence value of 0.53 Mq. In fact, with 
the Origlia et al. (2014) formula for mass loss the TP- 
AGB star is predicted to experience 9 thermal pulses be¬ 
fore leaving the AGB, while in all other COLIBRI mod¬ 
els the total number of thermal pulses is ~ 3 — 4. The 
corresponding lifetime is therefore longer, ~ 2.1 Myr, 
compared to ~ 0.5 — 0.9 Myr for the set of TP-AGB 
calcul ations made with the modified ISchrdder fc Guntzl 
(|2005|) relation. 

Even though the method to derive the mass-loss rates 
for RGB star s is int rinsically similar to that adopted by 
Origlia et al. (2014), i.e. dust radiative transfer calcula- 
tions that best tit the spectral infrared excess, the predic¬ 
tions for R GB mass loss obtaine d with the semi-empirical 
relation of Groenewegen (2012) are completely different. 


3. THE EPOCH OF MASS LOSS 

Through an analysis of the distribution of magnitudes 
and positions of the young white dwarfs in the core of 
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TABLE 2 

Wind Parameters for the Various PARSEC/COLIBRI 
Models. 


»;rgb 

VuiSC 

A^hb 

\Me] 

AfwD 

[Me] 

A£t2T 

[Myr] 

AtcHB 

[Myr] 

iwD 

[Gyr] 

0.1 

0.7 

0.84 

0.54 

107 

94.5 

10.98 

0.2 

0.5 

0.79 

0.54 

108 

94.9 

10.98 

0.3 

0.3 

0.72 

0.53 

111 

97.9 

10.98 

0.4 

0.13 

0.65 

0.53 

109 

95.4 

10.98 

Orl4 

Orl4 

0.69 

0.57 

111 

97.9 

10.98 

Grl2 

0.8 

0.88 

0.54 

108 

96.0 

10.98 


References. Orl4 

Orielia et al. 

1201411 

Grl2 

Groenewegen 

p0l2 



47 Tucanae, |Heyl et al.| (|2015|) measured the rate of dif¬ 
fusion due to relaxation in the cluster. The right panel of 
Fig. [^depicts the radial distribution of the brighter and 
fainter white dwarfs. Both groups of white dwarfs have 
the nearly same mass because they formed from stars of 
similar mass. The mass of the turn-off in 47 Tucanae 
varies by about 0.3% over 100 Myr, and the mass of the 
resulting white dwarf would vary by 0.05% over this same 
period. However, the fainter ones were formed earlier, so 
they have been diffusing through the cluster for a longer 
time. We have also added the radial distribution of the 
stars on the upper main sequence. This distribution is 
only slightly more concentrated than the young white 
dwarfs, indicating that there has been very little time for 
the young white dwarfs to have diffused through the clus¬ 
ter since their progenitors lost mass. Furthermore, the 
distribution of the upper-main-sequence stars is nearly 
identical to that of the red-giant branch stars even when 
we focus on seven-hundred stars near the tip of the red- 
giant branch. Usi ng the white-dwarf formation rate from 


Heyl et al. (20151 and the PARSEC models depicted in 
Fig: nr this corresponds to the last hundred million years 
along the red-giant branch, indicating that little mass 
loss has occurred up to about 100 Myr before the tip of 
the RGB. In particular one can pose the question how 
much diffusion has occurred between the upper main- 
sequence or equivalently the RGB and the young white 
dwarfs. 

We can play the diffusion back in time to what is pre¬ 
sumably the initial radial distribution of the stars, that 
of stars on the upper main sequence or the giants as 
depicted in Fig. [T] We assume a cooling curve for the 
white dwarfs and the best-fitting two-Gaussian model 
for the distribution of the white dwarfs as they diffuse; 
all that we vary is the time for main-sequence radial dis¬ 
tribution to diffuse and to appear like the distribution of 
white dwarfs younger than 50 Myr. We fit for the best¬ 
fitting timescales by performing a Kolmogorov-Smirnov 
test between the diffusion model and the distribution of 
stars. Fig. depicts the best-fitting diffusion models 
along with the observed distributions. If we assume a 
cooling curve for the white dwarfs as described in § [^ 
and a two-Gaussian model for the diffusion, we find that 
the mass-loss event coincides with the tip of the AGB, 
and the mass loss greater than 0.2 solar masses earlier 
than 20 Myr before the tip of the AGB can be excluded 
with ninety-percent confidence. A major mass-loss event 
on the RGB (greater than 0.2Mq can be excluded at the 
greater than the four-sigma level {p < 10“®) from the 
diffusion of the youngest white dwarfs alone. 



Fig. 2.— Radial distribution of the upper-main sequence stars 
as depicted in Fig. ^ and the white dwarfs younger than 50 Myr. 
There are 367 young white dwarfs, and their median age is 18.9 Myr 
since the tip of the AGB. The diffusion time between the two model 
distributions is 21.9 Myr. 


Fig .[^and|^also depict the Kolmogorov-Smirnov prob¬ 
ability obtained by evaluating the diffusion time-scale 
required to go from the best-fitting model for the UMS 
stars to the best-fitting model for a sample of young WDs 
(median age of 18.9 Myr since the tip of the AGB) as de¬ 
picted in Fig. This brings the epoch of major mass 
loss practically coincident with the TP-AGB. This prob¬ 
ability assumes that the theoretical diffusion model is 
fixed; in particular, it does not include the approximately 
ten-percent uncertainty in the diffusion timescale. This 
would shift the peak of the probability about 2 Myr in 
either direction. The statistical standard deviation in de¬ 
termining the median age of the white-dwarf sample is 
1.2 Myr. The theoretical cooling curve itself also yields 
an additional uncertainty in the diffusion timescale of 
about 4 Myr. This is obtained by comparing of the the- 
oretical cooling curve w ith the empirical cooling curve 


(Goldsbury et al. 

2012 

). Gombining these uncertainties 

m quadrature yiel 

ds the result that the mass-loss should 


have taken place 3.0±8.1 Myr before the tip of the AGB. 

The indications from the diffusion modelling show that 
the bulk of the mass loss is likely to have occurred while 
the white dwarf progenitors ascended the asymptotic gi¬ 
ant branch. The distribution of the young white dwarfs 
closely resembles that of the upper main sequence stars 
because they have not had time to relax to their new 
mass. It is unclear whether all of the mass loss occurs on 
the asymptotic giant branch or just the bulk of it. From 
these dynamical arguments it is likely that the mass-loss 
evolution calculated from the MESA model with equal 
values of rj depicted by the lowermost curve in Fig. 
is not correct. In this case during the first mass loss 
event, the mass of the progenitor decreases by 0.20 so¬ 
lar masses on the red-giant branch and 0.17 solar masses 
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Fig. 3. — Mass of a star that will become a modern-day white 
dwarf in 47 Tucanae as a function of time before its peak luminosity 
for the MESA models. The first epoch of mass loss in the model is 
as a red giant and the second is as an asymptotic giant star. The 
lowermost curve traces the mass loss for = I7AGB = 0.46 and 

the upper curves trace the mass loss for ?7rgb = 0.1, 0.2 and 0.3 
from top to bottom with the other quantities given in Table[^ The 
magenta curve depicts the KS probability of the epoch of mass loss 
obtaining calculating the diffusion of the UMS stars over a given 
time interval and comparing the resulting radial distribution with 
that of the young white dwarfs (median age 18.9 Myr) — the best 
fitting time for the mass loss to have occurred coincides with the 
TP-AGB. 


on the asymptotic giant branch. Similar considerations 
apply to the PARSEC-COLIBRI models as depicted in 
Fig. with ryRGB = 0.4. In this case the mass lost on the 
red-giant branch, 0.25 solar masses, is even larger, by a 
factor of two, than the mass expelled on the asymptotic 
giant br anch, 0.12 sol a r mas ses. The model computed 
with the Origlia et al. (2014) formalism predicts similar 
amounts of ejected masses, with the notable difference 
that the hnal mass and the TP-AGB lifetime are larger. 
Two of the diffusion models point to the AGB as the 
location of the mass loss, so it is likely that even more 
than two-thirds of the mass loss occurs then. 

The conclusion is that in 47 Tuc the mass of the 
horizontal-branch should be not much differ ent from that 
the main-sequence stars. Miglio et al. (2012) argued from 
asteroseismic observations of red giant and red clump 
stars in NGG 6791 that at least in this metal-rich cluster 
little mass is lost on the red-giant branch, less than a 
tenth of a so lar mass. If we use the value of rj = 0.2 as 
suggested by Miglio et al. (2012) the mass of the star de¬ 
creases more on the AGB as traced by the second curve 
from the top in Fig.[^ Here the star only loses 0.08 solar 
masses on the red-giant branch and 0.28 solar masses on 
the AGB in better accordance with the diffusion mod¬ 
els. The predictions of mass loss with the PARSEG- 
GOLIBRI model are of 0.11 and 0.25 solar masses, re¬ 
spectively for 77RGB = 0.2 and rymSC = 0.5. Further sup¬ 
port to the findings of our diffusio n models comes from 
the RGB mass loss rates derived by Groenewegen (20121 


Fig. 4.— Mass of a star that will become a modern-day white 
dwarf in 47 Tucanae as a function of time before its peak luminos¬ 
ity for the COLIBRI models. From top to bottom the four curves 
(solid lines) trace the mass loss for i?RGm = 0.1, 0.2, 0.3, and 0.4 
with the other qu antities given in Tab le The short-dashed line 
corresponds to the|Origlia et al.|||2014|| formulations for mass loss, 
while the up permost long-dashed line refers to the mass-loss rates 
pred icted bvIGroenewegenI||20 12|| for the RGB and with the mod¬ 
ified [Sdirdde^^^^uSzfj^DOSfTelat ion for the AGB. Note that in 
both cases the RGB rnass loss is derived from the results of the 
semi-empirical infrared-excess method, with quite different results. 


from the measured infrared excess in a sample of field 
RGB stars. In fact, adopting his proposed formulation 
we obtain that only « 0.02 Mq is expelled on the RGB, 
so that the star reaches the AGB with a mass that is in 
practice the same it had on the main sequence. In this 
case « 0.34 Mq is lost on the AGB. We emp hasize that 
the agreem ent was obtained with the original |Groenewe-| 


gen (2012) relation for the RGB mass loss rates, without 
introducin g any efficiency para meter. Moreover, as dis¬ 
cussed by Groenewegen (2012), the mass-loss rates de¬ 


rived in his study are consistent with the chromospheric 
estimates for RGB stars. It also appears that, surpris¬ 
ingly, the winds in metal-rich stars may operate similarly 
to those in more metal-poor stars such as in 47 Tuc. 

The radial distribution of the horizontal-branch stars 
themselves may help to constrain their masses. The core¬ 
helium burning phase on the horizontal branch lasts for 
about 100 Myr and stars will linger in the region of GMD 
that we define as the horizontal branch for about 80 Myr; 
this is more than a relaxation time of 30 Myr (Heyll 


et al. 2015), so if there is a significant amount of mass 
loss before the horizontal branch, one would expect the 
horizontal-branch stars to be less concentrated than the 
upper main sequence. Fig. shows that the horizontal 
branch of the evolved blue-straggler stars could contam¬ 
inate the HB region of the GMD due to saturation. To 
remove the saturated stars from the sample we have im¬ 
posed a mild cut on the uncertainty in the magnitude es¬ 
timates. Fig. also depicts the radial distribution of the 
giants and horizontal-branch stars. The giant stars have 
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a similar distribution to the upper-main sequence, and so 
do the horizontal-branch stars. If the horizontal-branch 
stars and the upper-main sequence stars were drawn from 
the same distribution, one would expect to find a devia¬ 
tion in the cumulative distribution as large as observed 
more than one third of the time (using the Kolmogorov- 
Smirnov test). Furthermore the number of horizontal 
branch stars after correcting for incompleteness (563) is 
in loose agreement with the theoretical duration of the 
horizontal branch of 80 Myr and the birth rate of whit e 


20151. 


dwarfs in the field of about 7 Myr ^ (|Heyl et al. 

Fig. [T] also depicts the radial distribution ol main- 
sequence stars whose mass is about O.65M0 as deter¬ 
mined from the PARSEC isochrone (in lavender); the 
distribution nearly coincides with the faint white dwarfs 
(in green). If the mass of the horizontal branch were 
C).65M0, this is the radial distribution to which the 
HB branch stars would ultimately relax. Of course, 
this would take about 100 Myr to reach completion, 
and the average age of a horizontal branch star is only 
40 Myr but the bulk of the relaxation would occur within 
40 Myr. The null hypothesis that the horizontal branch 
and the main-sequence stars of a mass of about O.65M0 
are drawn from the same radial distribution can be re¬ 
jected by a Kolmogorov-Smirnov test at the six-sigma 
level (p « 3 X 10“^^). 

By comparing the distribution of upper-main-sequence 
stars, giant stars, horizontal branch stars, young and 
old white dwarfs, we find that it is most likely that the 
stars that are currently evolving to form white dwarfs 
lose more mass as asymptotic giant branch stars than 
as red-giant-branch stars. In the context of the Reimers 
and Blocker or modified Schroder & Cuntz models for 
wind mass loss, parameters such as prgb ~ 0.1 and 
?7AGB ~ 0.7 best account for the observed radial dis¬ 
tributions of the stars. Conversely, all wind descriptions 
that predict the bulk of mass loss on the RGB conflict 
with the indications presented in this work. 

4. CONCLUSIONS 

We must emphasize that the results presented here 
proceed from two independent observations. The first 
is an inference of the time-scale of the major mass loss 
which rules out the RGB in favor of the AGB phase. 
This conclusion is nearly independent of the assumed 
distance to 47 Tucanae and the models of the horizontal- 
branch evolution. The time scale that we derive rests 
on three independent arguments: the dynamical relax¬ 
ation time from theoretical considerations, the white- 
dwarf cooling timescale and the duration of the red- 
giant-branch evolution that can be used to estimate the 
ages of the whit e dwarfs without reliance on white-dwarf 
cooling models (Beyl et al.||2015 1 — we count the num¬ 
ber of red-giant stars m the CiVID and determine the 
white-dwarf birthrate from the theoretical evolu tionary 
timescale throug h this portion of the CMD (see Golds- 
bury et al.|2012 for further details). All three agree and 
point to the AGB as the origin of the mass loss. 

Furthermore, even a large relative error in these 
timescales would only result in a slight change in the 
timescale of the mass loss relative to the duration of 
the horizontal branch because we are measuring the time 
elapsed between the mass-loss event and the appearance 
of a young ~ 20 Myr white dwarf. We find this to be 


about 20 Myr, so our estimate of the diffusion timescale 
would have to be underestimated by a factor of five to 
place the bulk of the mass loss on the red-giant branch. 
In this case it would be difficult to account for the diffu¬ 
sion between the bright and faint white dwarfs unless the 
white-dwarf cooling timescale were also underestimated 
by a factor of five, and the white-dwarf birthrate were 
overestimated by a factor of five as well. Because the 
white dwarfs are produced through stellar evolution, a 
revision of the white-dwarf birth rate would require a re¬ 
vision of the timescales for the entire post-main-sequence 
stellar evolution to achieve consistency with the numbers 
of stars observed in the CMD of 47 Tuc. In particular 
this would also increase the duration of the horizontal 
branch by about a factor of five as well, leading to fur¬ 
ther theoretical and observational inconsistencies. 

The second line of evidence rests on the observed ra¬ 
dial distribution of horizontal branch stars that strongly 
resembles the distribution of the upper-main sequence 
stars. The horizontal branch lasts long enough to suffer 
from diffusion if a significant amount of mass loss oc¬ 
curs before it. By using a mild cut in magnitude error 
we have eliminated the most saturated stars from our 
horizontal-branch sample. As shown in Fig.j^ these sat¬ 
urated stars are most likely the descendents of the blue- 
straggler stars. The number of horizontal branch stars 
in the field is also in accord with theoretical models and 
rate of formation of white dwarfs in the field. 

Our conclusions face two possible diffi culties: the pres¬ 
ence of multiple populatio ns in 47 Tuc (Anderson et al. 


2009 Milone et al. 


2012 b) and the conclusion from the 


current state of the art of horizontal-branch modelling 
that more mass loss is req uired to account for the ob - 
served horizontal branch (Idi Criscienzo et al. 2010). 


Milone et al. (2012b) found that the second geiieration 
dominates the population most strongly in the center of 
the cluster and the ratio of the two populations is con¬ 
stant with radius within the error bars where our ob¬ 


servations focus. Although Richer et al. (2013) found 
dynamical signatures and a radial gradient m the two 
populations, those conclusions were based on observa¬ 
tions far from the cluster core. In that outer field the 
relaxation time is much longer than in the core, so these 
initial differences have not yet been erased. The short 
dynamical time in the core compared to the age of the 
cluster ensures that dynamically the two populations be¬ 
have similarly; furthermore, the contribution of the first 
population is most modest in the core. 

On the sec ond front, recent synthetic h orizontal branch 
models (e.g. di Criscienzo et al. 2010) argue that the 
mass of the horizontal branch stars is significantly less 
than that of the mass sequence stars by about O.27M0. 
Bowever, a more firm conclusion on theoretical grounds 
may come only considering the simultaneous matching of 
the whole CMD of the cluster and exploring the possible 
degeneracy between various parameters (i.e. helium con¬ 
tent, metal mixture, and mass-loss efficiency). A detailed 
comparison both along the main sequence and through¬ 
out the post-main sequence evolution (especially the BB 
and AGB) with the data in all available bands including 
the ultraviolet would test the models and either support 
or contradict the conclusions drawn here from dynamical 
evidence. 

A key piece of evidence that could bolster these argu- 




























ments would be a direct measurement of the masses of 
the red giant and horizontal-branch stars in 47 Tucunae , 
perhaps through asteroseismology as Miglio et al. (20121 


did in NGC 6791 or by radial velocity measurements of 
binaries that include stars in these evolutionary stages. 
In any case the findings of this paper on 47 Tucanae are in 
clear contrast with other independent studie s that indi¬ 
cate the RGB as the phase of major mass loss (McDonald 
& Zijlstra 2015 Origlia et al. 2014). At the same time, 
they will also set a strong challenge to stellar evolution 
models, especially with regard to the detailed reproduc¬ 
tion of the morphology of the CMD (in particular the 
HB) of this cluster which is known to host two stellar 
populations with peculiar che mical mixtures and s lightly 
different helium abundances (Milone et al.|[^12a|. 
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